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ABSTRACT 

An inside-out model for the formation of haloes in a hierarchical clustering scenario 
is studied. The method combines the picture of the spherical infall model and a mod- 
ification of the extended Press-Schechter theory. The mass accretion rate of a halo is 
defined to be the rate of its mass increase due to minor mergers. The accreted mass is 
deposited at the outer shells without changing the density profile of the halo inside its 
current virial radius. We applied the method to a flat ACDM Universe. The resulting 
density profiles are compared to analytical models proposed in the literature, and a 
very good agreement is found. A trend is found of the inner density profile becoming 
steeper for larger halo mass, that also results from recent N-body simulations. Addi- 
tionally, present-day concentrations as well as their time evolution are derived and it 
is shown that they reproduce the results of large cosmological N-body simulations. 
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1 INTRODUCTION 

Numer ical s tudies (Quinn, Salmon fc Zur ek lll986t) : Frenk 
et al. Jl988f) ; Dub inski & Galberg <199ll) : C rone, Evrard 
& Richstone (Il994l) ; Navarr o, Fre nk & White (Il997l) . here- 
after NFW; Cole & Lacey Jl996ft; H uss. Jain & Ste inmet z 
119991) : Fukushige & Makino J 19971): Moore et al. Jl99d) . 
herea fter M GQSL; Jing & Suto feOOOl) . here after J S: Hern- 
quist fo)9C|), hereafter H90, Kravtsov at al. <ll998f) . Klypin 
at al. j200 If)') show that the density profiles of dark matter 
haloes are fitted by models of the form 



Pf(r) = 



{r/r s ) x [l + {r/r s y] v 



(1) 



More specifically NFW proposed a model with A = = 
l,v = 2, H90 proposed A = 1, fj, = 1, u = 3, MGQSL A = 
1.5, n = 1.5, v — 1 while JS proposed A = 1.5, fi = 1, v = 1.5. 
The logarithmic slope 7, of the density profile, defined 

by 



7(7-) 



dlnp f (r) 
dr 



= A + \w 



1 + {r/rsY 



(2) 



is, in all the above models, a decreasing function of radius. 
It is smaller than 2 near the centre of the system and larger 
near its virial radius. 
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An interesting quantity that characterizes the shape of 
the density profile is the concentration. This is defined by 
c = R v ir/r 3 where R v i T is the virial radius of the system. In a 
hierarchical clustering model, the concentration is a decreas- 
ing function of the mass of the system. Thus, smaller sys- 
tems, that formed earlier, have higher concentrations than 
larger ones. This reflects the high density of the Universe at 
the epoch of their formation. The concentration of the re- 
sulting structures is studied in detai ls by l arge N-body simu- 
lations. In particular, Bullock et al. (|2001j), BKPD hereafter, 
constructed a toy-model that describes accurately the time 
evolution of the concentration in a way consistent with the 
results of their large cosmological simulations. 

Although numerical experiments are the most power- 
ful method to study the formation of structures, the de- 
velopment of analytical or semi-numerical methods is very 
important as well since they help to improve our under- 
standing about the physical processes during the formation. 
Density profiles of equilibrium cold dark mat ter ha loes are 
studied by such methods (e. g. Sye r & White il998l) . Avila- 
Reese, Firmani & Hernadez (I1998l ). Raig, Go nzalez -Casado 

Nusser 
2002)). 

Modifications of the exten ded Pres s-Sch echter theory 
(PS) (Press & Schecht er dl974l) : Bower Jl99ll) : Bond et al. 
(1991), Lacey & Cole (1993)) based on the distinction be- 
tween minor and major me rgers (Manrique & Salvador-Sole 
(Il996l) : Ki tayam a & Suto (Il99d) : Salva dor-S ole, Solanes & 
Manrique dl99Sl) . Cohn, Bagla & White ll200ll) ) are also used 
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for sty dying the formation of structures. Recently, such a 
modified PS approximation was combined with a sp herica l 
infall model picture of formation by Manrique et al. (2003), 
MRSSS hereafter. Their results are in good agreement with 
those of N-body simulations. 

In this paper we use the formalism of MRSSS, with 
justified modifications, and the same model parameters as 
in BKPD. We compare the characteristics of the resulting 
structures with those in N-body results. In Section 2, 
we discuss the modified PS theory and its application to 
the calculation of the density profile. In Section 3, the 
characteristics of the resulting dark matter haloes are 
presented. A discussion is given in Section 4. 



2 EXTENDED AND MODIFIED 
PRESS-SCHECHTER THEORY 

One of the major goals of the spherical infall model is the 
PS approximation. It states that the comoving density of 
haloes with mass in the range M, M + AM at time t is given 
by the relation: 



N(M, t)AM = 



2 S c (t) p b0 
Tva 2 (M) M' 



do-(M) 
AM 



AM (3) 



where a(M) is the present-day rms mass fluctuation on co- 
moving scale containing mass M and is related to the power 
spectrum P by the following relation 



\M) = — / W 2 (kR)P(k)k 2 dk 



(4) 



where W is the Fourier transform of the window function 
used to smooth the overdensity field. The mass M and the 
radius R are related by the equation 



I „ 3 fi TO 0-ffo r>3 



M = rPb0 R ^ 



-R° 



(5) 



where pto is the present-day value of the density of the un- 
perturbed Universe, fl m o is the present-day value of the den- 
sity parameter (defined at any scale factor a by the relation 
Qm(a) — 8iTGpb(a)/(3H 2 (a))) and Ho is the present-day 
value of Hubble's constant, H. 

The only time dependent term of Eq. 3 is 8 c (t) that 
is the linear extrapolation up to the present epoch of the 
primordial density that collapses at t. It is calculated using 
the following arguments: In a model universe with cosmo- 
logical constant A, the radius r of a sphere having initial 
overdensity Ai, evolves according to the equation 



(6) 



where s = r/r» and r% is the initial radius. Hi is the value 
of the Hubble's constant at the initial time ti and g is given 
by the relation 

g(s) = Q mii (l + A l )(s- 1 - 1) + n A ,,( s 2 - 1) + 1 - |/iAi(7) 

and SIa,i is the initial values of the quantity 17a (a) = 
A/(3H 2 (a)). Eq.7 is derived under the assumption that the 
initial velocity Vi of the shell is w, = HiTi — v peCt i where 
the initial peculiar velocity, v pec ,i, is given according to the 



linear theory by the relation v peCy i = ^HiVifiAi (Pee- 
bles |l980)). A very good approximation of the fac tor fj is 
/i = fi m , i - 6 + ^[l-fim,i(l + O m ,i)] (Lahavetal. l|l991^ 1. 
The radius of the maximum expansion is rt a = st a ri, where 
sta is the root of the equation g(s) = that corresponds 
to zero velocity (As /At = 0). The sphere reaches its turn- 
around radius at time 



Ua ~ Hi 



i(s)ds 



(8) 



and then collapses at time t c = 2tt a - 

The scale factor a of the Universe obeys the equation 

da 

d7 

where 



H X?(a) 



(9) 



X(a) = 1 + fi m , (a 1 - 1) + n A ,o(c 



1) 



(10) 



and the subscript denotes the present-day values of the 
parameters. However, the time and the scale factor a are 
related by the equation 



t 



1 

Ho 



X-i(u)Au 



(11) 



Setting t — t c in the above equation and solving for a one 
finds the scale factor a c at the epoch of collapse. If we call 
Si, c (t) the initial overdensity required for the spherical region 
to collapse at that time t and take into account the linear 
theory for the evolution of the matter density contrast 5 = 
5p/p, we have 



5 cx 



1 X 1 / 2 (a) 
^ a 



X 



-3/2 



(u)Au = D(a), 



(Peebles 1980), then S c (t) is given 

*«(*) = M*)^j, 

where to denotes the present epoch. 
Usually 5 C is written in the form 



5 c (t) = 5i, c (t) 



D(t) D(t ) 
D(U) D(t) 



= <5crit(t) 



D(to) 
D(t) 



(12) 



(13) 



(14) 



where 5 C rit(£) is the linear extrapolation of the initial over- 
density up to the time t of its collapse. In an Einstein-de 
Sitter universe (fi m = l,f2 A = 0) this value is independent 
on the time of collapse and is <5 cr it ~ 1.686. In other cos- 
mologies it has a weak d epend ence on the time of collapse 
(e.g. Eke, Cole & Frenk ( 1996)). In a flat universe it can be 
approximated by the formula <5 crit (t) ~ 1.686f^°S 55 (t). 

The PS mass function agrees relatively well with the 
results of N -body simulations ( e.g. E fstathiou, Frenk & 
White lll985l). Efstathiou & Re es lll988l) : White, Efstathiou 
fc Fre nk Jl993l) . Lacey fc Col e lll994l) : Gelb & Bertschinger 
Jl994h : Bond & Myers \l99tt) ) while it deviates in detail at 
both the high and low masses. Recent impr oveme nts (Sheth 
& Tormen dl999l): Sheth, Mo & Tormen (1200 ll) . see also 
Jenkins et al. B2001TI ~1 allow a better approximation involving 
some more parameters. The application of the above approx- 
imation to the model studied in this paper is a subject of 
future research. 

Lacey & Cole ( |l993l) extended the PS theory using the 
idea of a Brownian random walk, and were able to calcu- 
late analytically tractable expressions for the mass function, 
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merger rates, and other properties. They show that the in- 
stantaneous transition rate at t from haloes with mass M to 
haloes with mass between M , M' + AM is given by 



r(M M',t)AM' = 



(I) 



2\V 2 dS c (t) 
n 



da(M') 



dt 



x exp 



a 2 (M') 
5l{t) 



dM> 



a 2 {M') 



a 2 (M) 

i 1 

a 2 (M') ~ a 2 (M) 



-i -3/2 



AM' 



(15) 



This provides the fraction of the total number of haloes 
with mass M at t, which give rise per unit time to haloes 
with mass in the range M' , M' + AM' through instantaneous 
mergers of any amplitude. 

An interesting modification of the extended PS theory is 
the distinction b etween minor and major m ergers (Manriquc 
& Salvador -Sole Jl996t) : Kitayama & Suto (Il996h: Sa lvador- 
Sole et al. Jl99St) : Percival, Miller & Peacock feOOch . Cohn 
et al. feOOlR 

Mergers that produce a fractional increase below a given 
threshold A m are regarded as minor. This kind of mergers 
corresponds to an accretion. Consequently, the rate at which 
haloes increase their mass due to minor mergers is the in- 
stantaneous mass accretion rate and is given by the relation 

pM(l + Am) 



r*(Af,t) = 



(M' - M)r(M -> M', t)AM' (16) 



Thus the rate of the increase of halo's mass due to the ac- 
cretion is 
AM(t) 



At 



C[M(t),t] 



(17) 



Before proceeding further with the model, it is useful to dis- 
cus briefly the cosmological considerations about the virial 
radius of a spherical system. Let A v i r (a) be the ratio of the 
overdensity of a sphere, that has collapsed and virialized 
at scale factor a, to the background density. This can be 
expressed by the form: 

p(a) 1 / " s 3 



A vir (a) = 



Pi, (a) 



— (-)( 



1 + Ai 



(18) 



where c/ is the collapse factor of the sphere defined as 
the ratio o f its fi nal radius to its turnaround radius. La- 
hav et al. lll99ll) applied the virial theorem to the viri- 
alized final sphere assuming a flat overdensity and found 
the collapse factor to be c/ w (1 — n/2)/(2 — n/2) where 
n — (Ar 3 a )/(3GM). For an Einstein-de Sitter Universe 
A v i r (o) « 187r 2 at any time. For flat models with cosmolog- 
ical constant, significantly goo d ana lytical approximations 
of Avir exist. Bryan & Norman (I1998T) proposed for A v i r the 
following approximation 



A vir (o) w (18tt - 82a; - 39x 2 )/n m (a) 



(19) 



where x = 1 — f2 m (a). 

MRSSS considered the following picture of the forma- 
tion of a halo: At time U an halo of virial mass Mi and virial 
radius Ri is formed and at later times it accretes mass ac- 
cording to the Eq.(17). Assuming that the accreted mass is 
deposited in an outer spherical shell without changing the 
density profile inside its current radius, then 

r-flvir(t) 



M vir (t) - Mi 



47rr p(r)dr 



(20) 



The current radius -R v ir contains a mass with mean den- 
sity A v i r (a) times the mean density of the Universe Pb(t). 
Therefore, 



Rvivit) 



3M vir (t) 



1/3 



(21) 



_47rA vir (t)p b (t) 
Differentiating with respect to t 

Mvir(t) A[ln[p b (t)A viY (t)]]~'~ 



p(t) = A vil (t)p b (t) 



r&.[M vll (t),t] 



At 



(22) 



Since one of the goals of this paper is the comparison of our 
results with the results of the N-body simulations of BKPD, 
we have used for A v i r (t) the same approximation as BKPD 
did -that is Eq.19- and not the constant value of 200 that 
MRSSS used. It is convenient to express Eq.22 in terms of 
the scale factor a instead of the time t. Thus, Eq.22 becomes: 



, . 3H 2 '(a)Q m (a) , . 
= o ^ A vir (a) x 



1 + 



M vir (a) 



8tyG 
3 dln[A vl r 
M v i T (a), a] \ a da 

where we used pb(a)a 3 = const, and 



(a)] 



(23) 



\M V 



AMa 
da 



iJ - 1 A- 1/2 (a)r, a n [M vir (f),t].(24) 



Integrating Eq.17 and using Eqs. 21 and 23, we obtain the 
growth of virial mass and virial radius and, in a parametric 
form, the density profile of haloes. 



3 DENSITY PROFILES OF DARK MATTER 
HALOES 

The results described in this section are derived for a flat 
universe with fl m ,o ~ 0.3 and JIa.o = 0.7. We used two 
forms of power spectrum. The first one -name d spec tl- is 
the one proposed by Efstathiou, Bond & White (1992). It is 
based on the results of the COBE DMR experiment and is 
given by the relation: 

Bk 

Ps P ectl(k) = ^ - ^ - (bfc)3/2 - (cfc)2] „ ]2/ „ (25) 

where a = (6.4/r)/i _1 Mpc, b = (3.0/r)ft _1 Mpc, c = 
(1.7/r)ft _1 Mpc and v = 1.13. Low-density Cold Dark Mat- 
ter (CDM) models in a spatially flat Universe (i.e. A > 0) 
are described for Y = flmfih. 

The second spect rum - named spect2- is the one pro- 
posed by Smith et al. <ll998T) and is given by: 

Ak n 

PsPeC " (fc) = [l + ai fcl/2 +a2fc + a 3 fc 3/2 +a4fc 2 ]6 ( 26 ) 

The values for the parameters are: n = 1, ai = 
-1.5598, a 2 = 47.986, a 3 = 117.77, a 4 = 321.92 and 
b = 1.8606. 

We used the top-hat window function that has a Fourier 
transform given by: 

3(sin(fc7?) - kRcos(kR)) 



W(kR) 



(kR) 3 



(27) 



The constants of proportionality A and B are found using 
the procedure of normalization for as = o~(R= 8h~ Mpc) = 
1. In Fig0 the resulting rms fluctuations for both spectra 
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log (M/M 




Figure 1. rms mass fluctuation as a function of mass. Solid line: 
for spectl, dotted line: for spect2 



Figure 2. For spect2: Density profiles for two haloes having 



the same mass 10 12 /i 
A m = 0.5 



Mq. Solid line: A T , 



0.21. Dotted line: 



are shown. It must be noted that we use a system of units 
with Munit = 10 12 /i _1 Mq, Runit = /i _1 Mpc and Unit = 



1.515 x 
1.5276. 



10 h years. In this system of units Ho/H un i t 



3.1 Present day structures 

In the approximation used in this paper, for given models of 
the Universe and power spectrum there is only one free pa- 
rameter, that is the value of the threshold A m (see Eq.16). 
We found that the resulting density profiles are sensitive 
to the value of A m . As an example, the density profiles of 
two systems with the same present-day mass 10 12 /i -1 Mq 
and different values of A m are plotted in Fig|2] Both den- 
sity profiles are derived from spect2. The solid line -shown 
in Fig|5J corresponds to the system derived for A m = 0.21 
while the dotted line to the system for A m = 0.5. The den- 
sity profile for smaller A m is steeper at both the inner and 
the outer regions. Additionally, for different A m the concen- 
trations of the haloes are different too (a detailed description 
of the way the concentration is calculated is given below). 
The system that results for A m = 0.21 has c v i r = 15.2, while 
the one for A m = 0.5 has c v i r = 8.7. In order to calculate the 
density profiles (that will be presented below), we used as a 
basic criterion the concentrations of the present-day struc- 
tures. In fact, we have chosen the values of A m = 0.23 and 
A m = 0.21 for spectl and spect2 respectively, because the 
concentrations resulting from these values are close to the 
results of the toy-model of BKPD. This model is constructed 
by BKPD to reproduce the results of their N-body simula- 
tions and it is able to give the concentration c v i r of a virial 
mass M v i r at any scale factor a. First, the scale factor a c 
at the epoch of collapse is calculated, solving the following 
equation 



a[M,(a c )] = a[FM vil (a)] 



(28) 



where F = 0.01 and M* is the typical collapsing mass. Then, 
the concentration is calculated using the formula 



Cvir.BKPD [M v i r (o) , O 



(29) 



where K — 4. We recall that the typical collapsing mass at 
scale factor a satisfies a[M,(a)] = 1.686D(1)/D(a). 

It is obvious that the above defined concentration de- 
pends only on the cosmology and the power spectrum used. 
Thus, for given cosmology and halo mass, the concentration 
Cvir, bkpd is known a priori without taking into consider- 
ation any particular form of halo growth. We applied this 
toy-model to find c v i r ,BKPD for the present-day structures. 

Another way to calculate the concentration is by using 
Cvir = Rvii/r2 where r-z is the radius where the logarithmic 
slope of the density profile equals 2. This radius is found by 
the following procedure: First, the resulting density profiles 
are fitted by the general formula of Eq.l. This is done by 
minimizing the sum 



S 



^2ip( r ) -Pfir) 



(30) 



where p c , r s , A, fi and v are fitting parameters. The min- 
imization is performed using the unconstrained subroutine 
ZXMWD of IMSL mathematical library. Then, r 2 is found 
by applying the following formula for n = 2 

\ 1 x l» 

t : r s (31) 

A + fj,u — n 



This formula gives the radius r n at which the logarithmic 
slope equals to n. According to the model presented in this 
paper, haloes grow inside-out. Thus, the value of c v i r repre- 
sents the way of halo growth. In Fig- 01 the concentration is 
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10° 10 1 10 2 

M . JM ., 

vir unit 

Figure 3. Concentration as a function of present-day virial mass. 
From the top, the first pair of curves are for spectl and the second 
for spect2. Solid curves: our results. Dotted curves: BKPD toy- 
model results. 



plotted as a function of the present-day virial mass. From 
the top of the figure, the first pair of curves (solid and dot- 
ted) correspond to spectf and the second pair to spect2. 
Solid curves show our results while dotted curves depict the 
results of the toy-model of BKPD. A very good agreement 
between the values of the concentration is shown. In partic- 
ular, concentrations resulting from spect2 are in agreement 
with those obtained for the model of BKPD for the whole 
range of mass presented. On the other hand, small differ- 
ences appear for very small and very large masses in the 
case of spectl. 

In Fig. |3] we present the density profiles of the re- 
sulting structures with present-day masses in the range of 
0.2 x lO 11 /i" 1 M to 8 x 1O 14 /i _1 M . The left-hand side fig- 
ures (al, bl, cl, dl, el) have been produced using spectl, 
while the right-hand ones using spect2. Figures (al) and 
(a2) correspond to mass 0.2 x 10 11 /i" 1 M Q , (bl) and (b2) 
have mass 10 12 /i _1 M Q , (cl) and (c2) to mass 10 13 h _1 MQ, 
(dl) and (d2) to mass 10 14 /i _1 Mq and (el) and (e2) cor- 
respond to mass 8 x 10 14 /i _1 Mq. Solid lines represent the 
resulting density profiles while dotted lines are the fits using 
the general formula of Eq.l. It is shown that the fits using 
the general formula of Eq.l are exact. We also fit every halo 
density profile using the analytical models that have been 
proposed in the literature (H90, NFW, MGQSL, JS) and are 
described in Section 1. The best fit of these models to our 
resulting profiles is shown in Fig. |3] (circles). This best fit is 
found by the minimizing procedure described above, for A, /i 
and v constants and equal to the proposed values, while p c 
and r B are the only fitting parameters. Best fit for the result- 
ing density profile in (al) is the H90 model, in (a2) and (b2) 
the NFW model, in (bl), (cl) and (c2) the MGQSL model 




-1.5 -0.5 -1.5 -0.5 



log(R/R vii .) log(R/R ir ) 

Figure 4. Density profiles as a function of radius. Solid curves: re- 
sulting density profiles. Dotted curves: fits of the resulting density 
profiles using the formula of Eq.l. Circles: best bit to our results 
using models proposed in the literature (H90, NFW, MGQSL, 
JS). Left-hand side: spectl. Right-hand side: spect2 



and in (dl), (d2), (el) and (e2) the JS model. Additionally, 
haloes of different mass are fitted well by different analytical 
models. This is due to the different inner and outer slopes 
of the density profiles. Inner slope, (defined as that at ra- 
dial distance r = 10~ 2 -R v ir), is an increasing function of the 
virial mass of the halo. For example, in the case of spect2 
the inner slope varies from 1.43 for M — 10 12 /i -1 Mq to 
1.65 for M = 8 x lO 14 ft _1 M . Additionally, outer slope -at 
r — R v ir- is a decreasing function of the virial mass and it 
varies from 3.67 to 2.64 for the above range of masses. 

Although density profiles resulting in simulations seem 
to be similar, systematic trends that relate them with the 
power spectrum have be en rep orted. For example, Subra- 
manian, Cen & Ostriker (2000) found in the results of their 
N-body simulations the following: for power spectra of the 
form P(k) oc k n the density profiles have steeper cores for 
larger n. Therefore, a dependence of the density profile on 
the power spectrum is expected. This dependence is shown 
in our results comparing the profiles of haloes with the same 
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800 



Figure 5. Exponent A as a function of present-day virial mass for 
both power spectra. It is shown that A is an increasing function 
of the virial mass. Solid curve: spectl. Dashed curve: spect2 



present day mass. It should be noted that the method stud- 
ied in this manuscript is applicable for the era of slow accre- 
tion when the infalling matter is in the form of small haloes 
that have mass less than Am times the mass of the parent 
halo. This kind of accretion occurs at the late stages of for- 
mation and thus determines the profile of the outer regions 
of the halo under study. However, the values of the inner 
slopes may be questionable. Real haloes have followed dif- 
ferent mass growth histories and thus their properties show 
a significant scatter about a mean value. Unfortunately, the 
method studied in this manuscript results one profile for a 
halo of given mass. Thus, its purpose is just to approximate 
the mean density profile of a large number of mass growth 
histories. Since the mass growth history resulting from the 
method is in good agreement with the mean growth his- 
tory resulting from N-body simulation -as it will be shown 
below- then the values of the inner slopes could be close to 
the ones of N-body simulations. A Monte Carlo analytic ap- 
proach based on the construction of a large number of mass 
accretion histories is under study. This study could answer 
to some of the above problems. 

In Fig. |S| the exponent A is plotted, that gives the 
asymptotic slope at R — > 0, derived by the general fit as 
a function of present-day virial mass for both power spec- 
tra. It is shown that the exponent A is an increasing function 
of virial mass. This trend of the inner density profile is also 
found in the results of recent N-body simulations (Ricotti 
(200l). 



lines show our resulting structures and the dotted lines show 
the mass growth curves of the model proposed by Wech- 
sler et al. |2002j). The curves of the left panel correspond 
to spectl while those of the right panel to spect2. From 
the top to the bottom, the curves correspond to masses 
2xlO 11 /i~ 1 M , lO la ft _1 M , lO 13 /i _1 M , 10 14 ^ -1 M Q and 
8 x 1O 14 /i -1 M respectively. It is obvious that massive 
haloes show substantial increase of their mass up to late 
times while the growth curves of less massive haloes tend 
to flatten out earlier. This behaviour of mass growth curves 
characterizes the hierarchical clustering scenario where small 
haloes are formed earlier than more massive ones. Addition- 
ally, it helps to define the te rm "fo rmation time" by a mea- 
surable way. Wechsler et al. define as formation scale 
factor a c the scale factor when the logarithmic slope of mass 
growth, (dlnM (a) / dlna) , falls below some specified value, 5. 
They use the value 5 — 2. It should be noted that this def- 
inition of formation scale factor differs from a c , defined by 
BKPD, since a c is the value of the scale factor at the epoch 
the typical collapsing mass is F times the virial mass of the 
halo. We found that the values of a c and a c for F = 0.01 and 
5 = 2 are different. This is also noticed in Wechsler et al. 
( 2002) since they state that S c and a c have similar values 
for 5 = 2 but for F = 0.015. However, the use of the value 
F = 0.015 in the toy-model of BKPD changes the resulting 
concentrations and so our basic criterion for the choice of 
the threshold A m is not satisfied. Therefore, it is preferable 
to choose a different value of 5 for the definition of a c , that 
of 5 = 1.5. In FigHJ the dotted lines show the ma ss gro wth 
curves of the model proposed by Wechsler et al. 120021) . In 
this model the mass growth is calculated using the relation: 



M vir (a) = M„ ir , exp[-a c 5(l/a - 1) 



(32) 



where M V i r< o is the present-day virial mass and the 
formation scale factor a c is defined by the condition 
dlnM(a)/dlna = 5 with 5 = 1.5. In Fig. |S| a very satis- 
factory agreement is shown, particularly for the less massive 
haloes. 

We have to note that our model haloes grow inside- 
out. Therefore, in early enough times -when the slope of the 
density is smaller that 2 all the way from the centre up to 
the current radius- it is meaningless to define c V i r . Once the 
building of the halo has proceeded beyond the point with 
slope 2, the evolution of c v i r is due to the growth of the 
virial radius and is given by 



r[M(a),a] = c vir (M 



R vil [M{a)} 



(33) 



where c v ; r (Mo) denotes the present-day concentration and 
R v i T [M(a)] and R v i T (Mo) are the values of the virial radius 
at scale factor a and at the present-day respectively. In Fig. 
[7|the time evolution of concentrations is plotted. Solid lines 
describe c v i r while dotted lines are c v ir,BKPD • More massive 
haloes have lower concentrations that evolve slower, while 
the concentrations of less massive haloes are higher and 
evolve more rapidly. 



3.2 Time evolution 

In Fig. |S| we plot mass growth curves. The curves show 
M v i r (a) as a function of a in a logarithmic slope. The solid 



4 CONCLUSIONS 

Since the formation of structures in a hierarchical clustering 
scenario is a complicated process, any attempt for the con- 
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log(a) log(a) 

Figure 6. Mass growth curves as a function of scale factor a. Left panel, right panel: spectl, spect2 respectively. From top to the bottom 
the curves correspond to masses 0.2 X 10 11 h~ 1 MQ, 10 12 /i _1 Mq, 10 13 /i _1 Mq, 1O 14 /i _1 M0 and 8 X lO 14 ft~ 1 M0 respectively. 




Figure 7. Concentration as a function of the scale factor a. Left panel: spectl. Right panel: spect2. Solid lines: our results. Dotted lines: 
the results of toy-model of BKPD. From top to the bottom the lines correspond to masses 0.2 X 10 1:l /i _1 Mq, 10 12 /i _1 Mq, 1O 13 /i _1 M0, 
10 14 h- 1 MQ and 8 X lO 14 h -1 M0 respectively. 



struction of analytical models requires a number of crucial 
assumptions. 

The model studied in this paper was proposed by 
MRSSS and assumes that 



(i) The rate of mass accretion is defined by the rate of 
minor mergers 

(ii) Haloes grow inside-out. The accreted mass is de- 
posited at the outer shells without changing the density pro- 
file of the halo inside its current virial radius 
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The first assumption indicates that structures presented in 
this paper formed by a gentle accretion of mass. The phys- 
ical process implied by the second assumption is that the 
infalling matter does not penetrate the current virial radius. 
This process requires an amount of non-radial motion. This 
amount has to be large enough so that the pericenter of 
the accreted mass is larger than the current virial radius. It 
should be noted that a density profile that results from a 
radial collapse has inner slope steeper than 2. It is the pres- 
ence of non-radial motion during the collap se tha t leads to 
inner slopes shallower than 2. (e.g. Nusser (2001), Hiotelis 
(2002), Subramanian, Cen & Ostriker (|200Ct lV Non-radial 
motions are always present in the structures formed in N- 
body simulations. 

Despite the above assumptions, the results of the model 
studied in this paper are in good agreement with the results 
of N-body simulations. The summary of these results is as 
follows: 

(i) Density profiles of haloes are close to the analytical 
models proposed in the literature as good fits to the results 
of N-body simulations. A trend of the inner slope of the 
density profile as an increasing function of the mass of the 
halo is also found, in agreement with recent results of N- 
body simulations. 

(ii) Concentration is a decreasing function of virial mass. 
Its values are in agreement with the results of numerical 
methods. 

(iii) Massive haloes increase their mass substantially up 
to late times. Growth curves of less massive haloes tend to 
flatten out earlier. The concentrations of less massive haloes 
evolve more rapidly while those of more massive haloes 
evolve slowly. 

Taking into account the number of assumptions and approx- 
imations used to build the model presented in this paper, we 
can conclude that the agreement with the results of N-body 
simulations is very good. Consequently, this model provides 
a very promising method to deal with the process of struc- 
ture formation. Further improvements to this model could 
help to understand better the physical picture during this 
process. 
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